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SECTION  I 


INTRODUCTION 

This  report  is  concerned  with  the  following  inverse  problem:  given  the 
strain  response  of  a  penetrator  vehicle  as  functions  of  time,  identify  the 
loads  as  functions  of  time.  The  objective  of  this  effort  is  to  use  experimental 
records  of  strains  in  vehicles  to  determine  the  forces  sustained  by  a  vehicle 
during  penetration.  In  this  effort,  attention  is  limited  to  the  development  of 
appropriate  techniques  for  the  problem  and  the  testing  of  certain  aspects  of  these 
techniques  by  comparison  with  solutions  and  experimental  results. 

The  identification  problem  for  off -normal  impact  is  quite  complex,  for  it 
involves  separate  flexural  and  axial  effects  that  arise  from  loads  distributed 
over  a  fairly  large  area.  In  the  present  investigation,  attention  is  restricted 
to  the  linear  response  of  the  vehicle,  that  is,  the  strains  are  assumed  to  remain 
small  and  the  response  is  assumed  to  be  elastic.  The  vehicle  is  modeled  as  a 
uniform  rod  with  the  nose  represented  by  a  rigid  body.  Contact  stresses  are 
represented  as  point  loads  and  moments  at  the  center  of  mass  of  this  rigid  body. 

The  identification  problem  is  subdivided  into  two  parts:  the  determination 
of  the  axial  force  and  the  determination  of  the  flexural  forces.  A  simple 
procedure  is  developed  so  that  if  the  axial  strain  is  known  as  a  function  of 
time  at  any  point  along  the  neutral  axis  of  the  vehicle,  the  axial  force  on  the 
nose  can  be  determined.  The  determination  of  the  flexural  forces  requires  that 
strains  be  measured  at  a  minimum  of  two  stations  along  the  length  of  the  vehicle; 
at  each  station,  strains  must  be  measured  at  three  points.  The  procedures  are 
outlined  in  Sections  2  through  4.  In  Section  5  some  sample  results 
are  given  which  illustrate  the  effectiveness  of  the  methods  for  the  axial  and 
flexural  problem.  The  methods  are  then  used  to  identify  the  forces  for  a  reverse 
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ballistics  experiment  (ref.  1).  The  results  are  evaluated  by  applying  these 
identified  loads  on  a  detailed  finite  element  model  and  comparing  computed 
strains  and  accelerations  to  experimental  records. 
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SECTION  II 


THE  FORCE  IDENTIFICATION  PROBLEM 


In  its  most  general  terms,  the  identification  problem  consists  of  deter¬ 
mining  the  distribution  of  normal  and  tangential  loads  on  the  vehicle  as  a 
function  of  time.  This  identification  must  be  accomplished  from  strain  gauge 
records  at  selected  points  in  the  vehicle  and  perhaps  one  or  two  accelerometer 
records.  The  elements  of  this  problem  are  illustrated  in  Fig.  1. 

When  considered  in  this  degree  of  generality,  the  identification  problem  is 
horrendously  difficult.  The  reason  for  this  is  explained  in  the  following. 

Consider  the  pressures  at  points  A  and  B,  shown  in  Fig.  1,  which  are  called  p^ 

B 

and  p^  ,  respectively.  The  response  of  a  strain  gauge  at  C  due  to  the  pressures 
A  B 

p  and  p  will  differ  because  the  signal  from  B  will  arrive  sooner  than  that 
n  n 

from  A.  However,  subsequent  to  the  arrival  of  the  wave  from  B,  there  is  no  way 
to  distinguish  whether  the  signal  arrived  from  A  or  B,  unless  a  strain  gauge 
is  located  between  A  and  B. 

It  is  possible  to  distinguish  the  signals  if  the  expected  depth  of  penet¬ 
ration  is  known  as  function  of  time.  However,  if  the  load  is  to  be  determined 
after  partial  penetration  of  the  nose,  this  information  is  of  little  help  since 
subsequent  to  penetration  beyond  B,  both  points  A  and  B  can  be  sources.  This, 
in  conjunction  with  practical  limitations  on  the  number  of  available  strain  gauge 
records,  renders  the  general  problem  of  distributed  load  identification  impractical. 

For  purposes  of  formulating  a  practical  force  identification  problem,  we 
will  limit  the  endeavor  to  net  resultant  forces  on  a  specified  segment  of  the 
vehicle.  The  following  assumptions  are  made: 


7 


Figure  1.  Depiction  of  general  force  identification  problem  for  a  penetrator 


^  m  HI  UlPyff,  (  I  I,  ^  I  t\.,+  >WM,«K«  <*«*  I(1M*  pinwMfcV  ,  (fl.^  !•«<#•«*  I1  A*»  II  «  H*  «*J  t!  *W*t»iflii*'4 


i.  Contact  between  the  vehicle  and  medium  occurs  primarily  in  the  nose 

region;  behind  the  nose,  the  medium  tends  to  separate  from  the  vehicle. 


ii.  No  significant  torques  about  the  vehicles  axis  are  exerted,  nor  does 
the  vehicle  spin  about  its  axis. 

The  first  assumption  is  based  largely  on  normal  impact  calculations  per¬ 
formed  by  CRT  (ref.  2) .  In  off-normal  impact,  the  DNA-SANDIA  (ref.  1)  experi¬ 
ments  have  shown  significant  backslap,  so  this  assumption  limits  the  method  to 
identification  of  early  loads.  Without  this  assumption  the  identification  problem 
would  become  very  complex.  The  second  assumption  is  based  on  the  observation 
that  in  impact,  no  mechanism  exists  which  would  exert  a  torque  on  the  vehicle. 


The  vehicle  is  modelled  as  a  uniform  elastic  bar  and  the  nose  represented 
by  a  rigid  body  of  finite  dimensions.  This  simplification  was  used  in  predicting 
the  response  to  normal  impact  of  the  DNA-SANDIA  vehicle  (ref.  3)  and  was  quite 
successful.  The  model  is  shown  in  Pig.  1. 

The  loads  to  be  identified  then  consist  of  the  following  forces  and  moments 
on  the  nose  (rigid-body  in  the  model) : 


i.  the  axial  force,  f  (t) ,  and  transverse  forces  f  (t)  and  f  (t) 

X  V  z 

ii.  the  moments  m  (t)  and  m  (t) . 

y  z 

In  order  to  describe  the  equations  governing  this  model,  it  is  convenient 
to  introduce  two  coordinate  systems; 

i.  a  moving  coordinate  system  (x,y,z)  embedded  in  the  projectile  such  that 
x  always  lies  along  the  axis  of  the  vehicle  as  shown  in  Fig.  2,  with 
the  origin  located  at  the  interface  between  the  rigid  body  and  the 
deformable  bar; 

ii.  a  coordinate  system  (X,Y,Z)  fixed  to  the  target,  which  is  considered 
inertial. 
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The  displacement  of  a  point  in  the  vehicle,  r,  relative  to  the  vehicle's 

coordinates  is  given  by  a  vector  v  with  components  v^,  v  ,  vz,  which  are  assumed 

to  be  small.  The  displacement  of  the  center  of  mass  of  the  vehicle  is  denoted 

-> 

by  V.  The  angular  velocity  and  angular  acceleration  of  the  vehicle  are  u)  and 
->■ 

a.  The  acceleration  of  any  point  of  the  vehicle  is  then 


••  *•  • 

R  =  u  +  2wxv  +  axr  +  (ox  (co  x  r) 

(2.1) 

u  =  V  +  v 


Since  we  assume  the  vehicle  is  not  spinning,  its  only  nonzero  angular 
velocities  are  (0^  and  <0^ .  Moreover,  the  vehicle  is  considered  slender  and  the 
rotatory  inertia  of  the  cross-section  will  be  neglected.  From  these  assumptions 
it  follows  that 

“>■  “  ->*  — 
r  =  x  i  x  =  x  -  xc 

(2.2) 

<o  =  (0  j  +  cok 

yJ  z 

The  acceleration  components  in  the  (x,y,z)  systems  are  given  by 


••  •  *  2  2  - 

R  =  vi+2(tov  -  co  v  )  -  (co  +  to  )x 
x  xvyz  zy  y  z 


R  =  u+2cov  +  ax 
y  y  z  x  z 


(2.3) 

(2.4) 


R  = 
z 


2(0  V 
y  x 


a  x 

y 


(2.5) 
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Linear  strength  of  materials  theory  will  be  employed  for  describing  the 
dynamic  deformation  response.  Simple  first  order  rod  theory  which  neglects 
the  effects  of  transverse  inertia  will  be  used  for  the  axial  response.  Euler 
beam  theory  will  be  used  for  the  bending  response.  The  cross-section  of  the 
vehicle  is  assumed  to  be  axisymmetric  so  that  the  moment  of  inertia  of  the  cross- 
section  is  isotropic.  Therefore  bending  in  the  x-y  planes  and  in  the  y-z  planes 
is  uncoupled  for  any  choice  of  the  orientation  of  the  y  and  z  coordinates. 
Moreover,  the  axial  deformation  is  uncoupled  from  bending  because  of  the  assump¬ 
tion  of  small  displacements. 

The  governing  equations  are  then 


2 

c  v  =  R 

X,XX  X 


(2.6) 


c 


c 


2 

b 


2 

b 


v  ,xxxx  =  R 

y  y 


V  =  R 

z,xxxx  z 


(2.7) 


(2.8) 


where  commas  denote  derivatives.  The  accelerations  are  given  by  Eqs.  (2.3)  to 
(2.5)  and 

p  =  density  of  vehicle  material 
E  =  Young's  modulus  of  vehicle  material 
A  =  cross-sectional  area 

I  =  moment  of  inertia  of  cross  section  about  both  the  y  and  z  axes 
2 

c  =  E/p,  elastic  rod  wavespeed 
c?  =  El/ pA 

D 


By  using  Eqs.  (2.3)  to  (2.5),  Eqs.  (2.6)  to  (2.8)  can  be  written  as 
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2 

c.v  -  u  =  b 

b  y,xxxx  y  y 


2  •• 
c,  v  -  u  =  b 

b  z,xxxx  z  z 


where 


•  •  '2.  2  — 
b  =  2(w  v  -  w  v  )  -  (w  +w)x 
x  y  z  z  y  y  z 


(2.10) 


(2.11) 


(2.12) 


b  =  2w  v 
y  z  x 


a  x 
z 


(2.13) 


b 

z 


-2a)  v  -ax 
y  x  z 


(2.14) 


The  terms  b  ,  b  ,  and  b  couple  Eqs.  (2.9)  to  (2.11).  In  order  to  uncouple 
x  y  z  ” 

the  equations  we  assume  that  these  terms  vanish.  Since  the  rigid  body  motion 
-> 

V  is  independent  of  x,y  and  z,  its  derivatives  vanish  and  we  can  replace 
Eqs.  (2.9)  to  (2.11)  by 


2 

c  u 

x,xx 


ii 

x 


0 


(2.15) 


2 

c,  u  -  u 

b  y,xxxx  y 

2 

c,  u  -  u 

b  z , xxxx  z 


0 


0 


(2.16) 


(2.17) 


The  displacements  u  and  u  and  their  second  derivatives  must  be  small,  but 
y  z 

u  may  be  arbitrarily  large,  although  its  first  derivative  must  also  be  small. 

X 

The  boundary  conditions  are 
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A  E  u  =  Mu  -  f  (t) 
X,X  X 


at  x  =  0 


A  u  =0 
x,x 


at  x  =  L 


-E  I  u  =  M  (ii  -  £vi  )  -  f  (t) 

y ,  xxx  y  y ,  x  y 


E  I  (u  -  £u  )  =  J  ii  -  m  (t) 

y,xx  y,xxx  y,x  z 


-E  I  u  =  M(u  -  )  -  f  (t) 

z ,  xxx  Z  Z  ,  X  z 


EI(u  -£u  )  =  jii  +m  (t) 

z,xx  z,xxx  zfx  y 


at  x  =  0 


u  —  u  =0 

z,xx  z,xxx 


at  x  =  L 


where  J  is  the  moment  of  inertia  of  the  rigid  body.  Initial  conditions  provide 
the  initial  velocity  of  the  projectile. 

The  strain  at  any  point  x,y,z  in  the  structure  is  given  by 


e(x,  y,  z,  t)  =  u  -  y  u  -  z  u 

x,x  J  y,xx  z,xx 


The  identification  problem  then  consists  of  the  following:  given  strain  histories 

at  selected  points  x.,  y.r  z.,  determine  the  forces  f  (t) ,  f  (t) ,  and  f  (t) 

i  -*1  i  x  y  z 

and  the  moments  ra  (t)  and  m  (t) .  For  convenience,  we  denote  by  (F  (t) }  the 

y  z 

complete  matrix  of  forces  to  be  determined:  {  f  ,  f  f  ,  m  ,  m  }. 

x  y,  z  y  z 


S'Jhen  strain  records  at  three  different  points  y^ ,  z^  are  given  at  any 
cross  section  x^,  then  Eq.  (2.21)  can  be  solved  uniquely  to  yield  the  response 
(r (t) }  =  {u  (x . , t) ,  u  (x . , t) ,  u  (x . , t) }  provided  that  the  matrix 

a  jl  y  / xx  i  z  f xx  1 


[A]  = 


1 

1 


■y 

■y 


(2.23) 


1 


-y 


3 


-2 


3 


is  not  singular 

The  matrix  [A]  is  singular  if  and  only  if  all  three  points  are  collinear. 

Thus  any  triangular  arrangement  of  the  points  y . ,  z .  is  sufficient  to  permit 

^  3 

the  determination  of  {r(t)}. 


Once  (r(t)}  is  given.  The  complete  identification  problem  can  be  subdivided 
into  three  uncoupled  parts: 

i.  given  e  =  u  (x.,  t) ,  determine  f  through  Eq.  (2.15)  and  boundary 

X|X  1  X 

conditions  (2.18) 

ii.  given  u^  ^  (x^,  t) ,  determine  and  fy  through  Eq.  (2.16)  and 
boundary  conditions  (2.19) 

iii.  given  u^  ^  (x^  t) ,  determine  my  and  fz  through  Eq.  (2.17)  and 
boundary  conditions  (2.20). 


Problem  (i)  is  an  axial  problem,  problems  (ii)  and  (iii)  identical 
flexural  problems;  they  are  discussed  in  Sections  3  and  4,  respectively. 

For  two  dimensional  problems,  where  the  motion  of  the  vehicle  is  planar,  only 
problems  (i)  and  (ii)  need  be  solved. 
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SECTION  III 


AXIAL  FORCE  IDENTIFICATION 


Vie  consider  here  the  problem  of  the  identification  of  the  axial  load.  The 
axial  response  problem  is  governed  by  the  equations  and  boundary  conditions 
(subscripts  x  are  omitted  in  this  section) 


U'xx  2  u 
c 


AE  u,  =  M  u 
x 


u,  =0 
X 


at  x  =0 


(3.1) 


(3.2) 


at  x  =  L 


which  correspond  to  Eqs.  (2.9)  and  (2.18).  The  initial  conditions  are  taken 


to  be 


u(x,0)  =  0 


u(x,0)  =  0 


(3.4) 


(3.5) 


The  problem  is  to  obtain  the  time  history  of  the  net  force  on  the  nose, 
f  (t) ,  given  the  axial  strain  time  history  at  one  point. 

We  will  use  Laplace  transform  techniques.  Let  the  Laplace  transform  of  a 
given  variable  be  denoted  by  a  corresponding  upper  case  letter,  i.e.  U(x,s)  is 
the  Laplace  transform  u(x,t),  so  that 


U(x,s)  =  f  u(x,t)e  St  dt 

J  o 


(3.6) 


. 


. . MU . Willi . . . I . Ill . Ill . Ill* . a . Ill . Bill® . I . I . .  ilHfttunii  ijHItitif . . . if . 


Similarly  E(s)  is  the  Laplace  transform  of  e(t),  F(s)  the  Laplace  transform  of 
f(t). 


Taking  the  Laplace  transform  of  Eq.  (3.1)  and  noting  that  the  initial 
conditions  are  zero,  g: 


U,  (x,s)  =  —  U(x,s) 

AA  4 

c 


(3.7) 


which  has  the  solution 


U(x,s)  =  + 


(3.8) 


The  Laplace  transforms  of  the  boundary  conditions,  (3.2)  and  (3.3),  give 


AEU,x(0,s)  =  M  s  u (0,s)  -  F(s) 


(3.9) 


U,x(L,s)  =  0 


(3.10) 


Substituting  Eq.  (3.8)  into  Eqs.  (3.9)  and  (3.10)  yields  two  linear,  algebraic 
equations  in  and  K2,  the  solution  of  which  is 


Kx  =  e2Ls/c  F (s)/G (s) 


(3.11) 


=  F (s)/G (s) 


(3.12) 


where 


G (s)  =  s  [(Ms  +  ^)  e2Ls/c  +  (Ms  -  ^)  ] 


(3.13) 
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From  Eq.  (3.8),  it  also  follows  that  the  Laplace  transform  of  the  strain  at 


any  point  x^  is  given  by 


£(Vs>  =  f  I-KJe-“l/c  + 


(3.14) 


Substituting  Eqs.  (3.11)  and  (3.12)  into  (3.14)  and  solving  for  F(s)  yields 


F(s)  =  iifer  ceui's)  G(s) 

H(s)  =  exis/c  -  e (2I-x±) s/c 

Substituting  Eq.  (3.13)  into  Eq.  (3.15), 
£(x. ,s) 

F(s)  =  —?■-  ( (Mcs  +  AE)e  /c  + 

H  (s) 


we  can  write  the  latter  as 


(Mcs  -  AE)) 


(3.15a) 


(3.15b) 


(3.16) 


The  time  domain  counterpart  of  Eq.  (3.16)  is 

x.  2L-x. 

f (t  +  —  )  -  f(t  +  - -}  =  Mce(x.,t)  -  AEe(x.,t) 

C  C  11 

+  Mce(x.  ,t  +  —  )  +  AEe(x.  t  +  —  )  (3.17) 

i  c  i  c 

Thus  by  shifting  in  time  by  (2L-x^)/c,  we  obtain  from  Eq.  (3.17)  that 

2L-x.  x. 

f(t)  =  f(t  -  2  (L  -  x.)/c)  +  AE  [£(x.,t  - -)  -  £(x  ,t  +  -±)  ] 

X  1C  1C 

x.  #  2L-x. 

-  Me  [e  (x.,t  +  —  )  -  e  (x.,t - -—)  ]  (3.18) 


From  the  above  equation,  it  is  apparent  that  once  an  axial  strain  record  is 
known  at  any  given  point,  the  force  history  can  be  determined.  Although  a 
Laplace  transform  approach  has  been  used,  from  the  time  shifts  that  occur  in 
the  equations,  it  is  apparent  that  this  result  follows  from  the  nature  of  the 
solution  to  the  wave  equation.  The  solution  always  consists  of  a  wave  induced 
by  the  force  which  moves  to  the  left  and  an  image  moving  to  the  right,  shifted 
in  time  by  L/c,  which  maintains  the  stress  free  condition  on  the  right  end  of 
the  rod.  The  Laplace  transform  approach  serves  to  formalize  these  results. 


SECTION  IV 


FLEXURAL  FORCE  IDENTIFICATION 


From  each  flexural  equation,  two  quantities,  a  force  and  a  moment,  must 
be  determined.  The  flexural  equation  and  associated  boundary  conditions,  Eqs. 
(2.10)  or  (2.11),  and  (2.19)  or  (2.20),  can  be  written 


1  .. 

W'xxxx  +  2  W 

Cb 


w,  =  c,  (w  -  £w,  )  -  f 
XXX  1  X 


at  x  =  0 


w,  -  Zw,  =  c_  w,  -  m  at  x  =  0 

xx  xxx  2  x 


w,  =0  at  x  =  L 

XX 


w,  =0  at  x  =  L 

XXX 


where 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


The  problem  is  to  identify  f(t)  and  m(t)  from  records  of  the  curvature 
wfxx  (x^,t).  As  in  Section  3,  we  wil?  use  a  Laplace  transform  approach  and 
denote  the  Laplace  transform  by  upper  case  letters. 


Taking  the  Laplace  transform  of  Eq.  (4.1)  yields 


W'x*»c+-2  W'° 

cb 


for  which  the  solution  is 


Y1X  -Y1X  Y2x  _Y2X 
W(x,s)  =  e  +  e  +  K3  e  +  K4  e 


Yi  -  »  +  s’VjF 


y2=  (1 


-  5>l/T 


and  K.  are  functions  of  s,  and 


i  =  /=1 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


The  boundary  conditions  at  x  =  0  and  x  =  L  can  be  written  in  partitioned  matrix 
form  in  terms  of  W(x,s)  as  follows 


] 


lH12] 


lH21]  lH22} 


K1 

F 

K_ 

M 

2 

► 

K 

0 

3 

K* 

0 

4J 

> 

(4.13) 


lHll3  = 


cls2d  -  Y1«  +  Yx3 
Y-,2  (1  -  Y-,£)  ~  c  s2Y, 


c^s2(l  +  yxW  -  Yj3 

2  2 
y,  (1  +  Yi  S-)  +  c_s  Y, 


(4.14a) 


If  the  curvature  is  given  at  two  points,  and  then  we  can  write 


W(x^,s) 


W(x2,s) 


[L11]  IL12] 


Prom  Eq.  (4-13),  it  follows  that 


,x  ...  .  N 


-IH22]  [H21]  \ 


Hence 


W(x^,s) 


W(x2,s) 


-  «12>  "W 


/ 

■x  w21)|  . 1 


If  we  call  the  above  square  matrix  [G] ,  Eqs.  (4.14)  and  (4.18)  yield 


-1  1  -1  W(x^,s) 

("ill  -  W12]  M22l  W21]  !g] 

/  W(x_,s) 


An  explicit  transformation  of  the  above  into  the  time  domain  is  impossible. 
Instead,  the  operations  of  Eq.  (4.19)  were  carried  out  in  the  frequency  domain, 
using  Fast  Fourier  Transform  algorithms. 

The  procedure  consists  of  the  following.  Transforms  are  obtained  of  the 

response  w,  at  two  points.  The  transfer  function  relating  the  input  to  the 
XX 

response  is  then  evaluated;  its  product  with  the  response  yields  the  transform 
of  the  forcing  function.  The  inverse  transform  of  the  forcing  functions  are 
then  taken  to  obtain  the  forcing  functions  in  the  time  domain. 

The  submatrix  [H^l  is  singular  at  s  =  0  (rigid  body  mode)  and  at  the 
natural  frequencies.  To  make  the  inversion  possible,  artificial  damping  is 
added  to  the  system.  This  displaces  the  zeroes  from  the  imaginary  axis  and  so 
avoids  the  singularities. 

Note  that  the  number  of  distinct  forces  that  can  be  identified  depends  on 
the  number  of  response  records.  Two  response  records  are  needed  to  identify  the 
force  and  moment  on  the  nose.  If  backslap  is  to  be  identified,  it  could  be 
represented  by  an  inhomogeneous  term  in  the  boundary  condition  Eq.  (4.5);  res¬ 
ponses  at  three  points  along  the  axis  of  the  vehicle  would  then  be  needed. 
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SECTION  V 


ILLUSTRATIVE  RESULTS 

The  force  identification  procedures  for  the  axial  force  and  the  flexural 
forces  have  been  tested  through  the  use  of  analytic  and  finite  element  solutions. 
Analytic  solutions  were  used  to  test  the  procedure  as  follows:  the  strain  was 
predicted  for  a  prescribed  forcing  function  at  the  nose,  and  this  strain  was 
then  used  as  input  to  the  identification  procedure.  The  identified  load  was 
then  compared  with  the  actual  force.  The  finite  element  test  was  similar  except 
that  finite  element  solutions  were  used  rather  than  an  analytical  solution  to 
obtain  the  strains  from  the  forces. 

In  both  tests,  the  vehicle  described  in  Fig.  3  was  used.  The  first  test 
involved  only  an  axial  load  with  the  time  history  shown  in  Fig.  4.  The  responses 
at  x  =  7in.  and  15in.  are  shown  in  Figs.  5  and  6.  Figure  7  shows  the  force  as 
computed  from  the  identification  procedure  given  in  Section  3,  and  compares  the 
results  to  the  applied  load.  As  can  be  seen,  the  identified  axial  load  corres¬ 
ponds  closely  with  the  actual  load.  There  are  some  discrepancies  at  the  points 
at  which  the  slope  of  the  actual  load  is  discontinuous;  the  identified  load  at 
those  points  tends  to  smooth  over  these  discontinuities. 

The  same  procedure  was  repeated  with  the  finite  element  solution.  The 
finite  element  solution  shows  the  effect  of  the  vehicle  idealization 
and  the  sensitivity  of  the  identification  procedure  to  spurious  oscillations. 
Finite  element  solutions  have  spectral  amplification  and  cut-offs  similar  to 
strain  gauges.  The  identified  force  is  compared  to  the  actual  force  in  Fig.  8. 

As  can  be  seen,  for  the  finite  element  solution,  the  identified  force  does  not 
correspond  as  closely  to  the  actual  input  force  as  in  the  case  of  the  analytic 
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solution.  There  is  a  more  pronounced  oscillation  about  the  actual  solution. 

Some  of  these  oscillations  may  be  removable  by  digital  filtering. 

In  the  third  example,  a  moment  with  a  time  history  was  applied  to  the  nose 
in  addition  to  an  axial  force  on  a  finite  element  model.  The  identified  loads 
are  compared  to  the  applied  loads  in  Fig.  9.  In  this  case,  the  identification  is 
not  as  accurate  as  for  the  axial  force.  The  use  of  numerical  transforms  introduces 
substantial  errors  unless  tremendous  resolution  is  used  in  the  frequency  domain. 

The  identification  shown  here  represents  a  compromise  between  adequate  resolution 
and  reasonable  computer  time.  The  maximum  moment  and  salient  aspects  of  the 
time  history  are  well  represented;  but  the  predicted  moment  exhibits  a  significant 
lag. 

Numerical  experiments  have  also  been  carried  out  in  identifying  both  the 
moment  and  shear  from  curvature  records  at  two  stations.  These  have  been  quite 
unsatisfactory.  To  obtain  a  transfer  function  which  is  invertible,  considerable 
damping  is  needed,  and  the  identified  forces  and  moments  are  quite  sensitive 
to  the  magnitude  of  the  damping.  Therefore,  in  the  following  section,  only  the 
transverse  force  was  identified  for  the  flexural  problem. 


TIME  (msec) 


Figure  8. 


Comparison  of  actual  and  identified  f (t)  using  strain 
at  x  =  7.0  from  finite  element  solution. 
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MOMENT 


240-1 


TIME  (msec) 


Figure  9.  Comparison  of  actual  and  identified  moment  using  strain 
from  finite  element  solution. 
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SECTION  VI 


LOADS  FOR  DNA  REVERSE  BALLISTIC 
EXPERIMENTS 

We  will  here  describe  the  loads  identified  for  the  DNA  reverse  ballistic 
experiments  conducted  at  Sandia  Laboratories  (ref.  1) .  The  study  is  limited  to 
Test  3  in  which  the  projectile  was  in  a  3  degree  nose-down  attitude;  i.e.  the  angle 
of  attack  was  3  degrees.  Test  data  consisted  of  four  accelerometer  records,  two 
axial  and  two  lateral,  and  seventeen  axial  strain  gauges. 

In  this  Section,  we  will  describe  the  simplifications  which  were  necessary 
to  generate  a  model  suitable  for  the  force  identification.  Once  the  forces  were 
identified,  they  were  applied  to  a  detailed  finite  element  model  of  the  projectile. 

The  strains  and  accelerations  predicted  by  the  finite  element  model  are  then 
compared  to  experimental  results. 

A  schematic  of  the  projectile  is  shown  in  Fig.  10.  The  idealization  employed 
for  the  force  identification  is  shown  in  Fig.  11.  Several  features  of  the 
projectile  had  to  be  omitted  because  the  identification  model  is  limited  to  a 
single  rigid  mass  and  a  rod  of  uniform  properties.  The  rigid  mass  was  used  to  model 
the  portion  of  the  projectile  between  0  and  3.45  inches.  The  uniform  rod  was  then 
given  the  dimensions  of  an  average  cross-section.  In  the  actual  model,  we  have 
four  distinct  sections  behind  the  nose  (1)  the  kennertium  plug  and  steel  barrel 
from  3.45  inches  to  7.45  inches;  (2)  the  1.702  inch  O.D.  section  from  7.45  inches 
to  9.64  that  houses  the  front  accelerometers;  (3)  the  tapered  section  from  9.64  inches 
to  13.04  inches,  in  which  the  O.D.  increases  from  1.702  inches  to  1.900  inches; 

(4)  the  last  constant  diameter  section  from  13.045  inches  to  18.15  inches.  The 
idealized  model  is  too  stiff  in  comparison  to  the  first  half  of  the  vehicle,  and  too 


flexible  in  comparison  to  the  second  half. 

The  loads  were  identified  from  the  strain  gauge  records  at  station  9.6  inches, 
gauges  3  and  10.  (shown  in  Figs.  17  and  18).  This  location  was  chosen  because  it 
is  as  far  as  possible  from  the  irregularities  near  the  front  accelerometer  and 
relatively  clean  records  were  obtained  in  Test  3.  The  force  was  identified  only 
for  0.4  msec;  after  that  time  the  vehicle  penetrates  beyond  the  nose. 

The  identified  axial  and  transverse  loads  are  shown  in  Fig.  12.  The  maximum 
axial  load  is  130,000  lbs.  and  occurs  at  300  microseconds;  the  maximum  transverse 
load  is  26,000  lbs.  and  occurs  at  50  microseconds.  Both  of  the  records  shown  were 
extensively  filtered. 

These  loads  were  then  applied  to  the  finite  element  model  shown  in  Fig.  11  for 
a  0.6  msec  simulation,  the  loads  at  0.4  msec  were  held  constant.  The  model  was 
solved  by  the  program  WHAMS  -  2D  (ref.  4)  using  explicit  time  integration.  All 
acceleration  results  were  smoothed  by  a  10  point  averaging  filter  (ref.  5).  The 
acceleration  results  are  compared  to  the  four  accelerometer  records  for  Test  3  in 
Figs.  13  to  16.  As  can  be  seen,  the  axial  accelerations.  Figs.  13  and  14,  agree 
quite  well  with  the  experiment,  indicating  that  the  axial  load  identification  was 
successful.  The  agreement  of  the  lateral  records.  Figs.  15  and  16,  is  quite  poor. 

Strain  gauge  output  is  compared  to  the  experiment  in  Figs.  17  to  23.  Particularly 
for  the  first  0.4  msec,  the  agreement  is  satisfactory.  The  increasing  discrepancy 
after  this  time  is  understandable  since  the  loads  were  only  identified  to  0.4  msec. 
However,  it  is  unclear  why  the  agreement  in  lateral  accelerations  is  so  poor  in 
comparison  to  the  strains.  A  possible  explanation  may  lie  in  the  nonuniqueness  of 
the  identification  problem.  If  a  different  model  were  selected,  the  nor.  vniqueness 
may  lead  to  substantially  different  lateral  forces  but  still  the  same  strains.  It 
is  this  nonuniqueness  which  is  also  undoubtedly  related  to  the  numerical  difficulties 
of  inverting  the  problem. 
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STA  7.45  ACC-Z 


EXPERIMENT 


STA  16.35  ACC-Z 


STA  7.45  ACC-R 


Figure  15.  Lateral  acceleration  (y-component)  at  station  7.45  (Accelerometer  1) 
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0.4  0.6 

TIME,  msec 


Figure  16.  Lateral  acceleration  (y-component)  at  station  16.35  (Accelerometer  3) 
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ATTN:  DRXSY-D,  J.  Sperrazza 

U.S.  Army  Missile  MD  Coosand 
ATTN:  SSiC 
ATTN:  F.  Fleming 
ATTN:  DSCPK-PE,  W.  Jann 

U.S.  Army  Mobility  Equip  RSD  Cmd 
ATTN:  DRK'E-KC 
ATTN:  ORDDKE-IS 
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DEPARTMENT  OF  IE  ARMY  (Continued) 

U.S.  Army  Nuclear  &  Chemical  Agency 
ATTN:-  Library 

U.S.  Army  War  College 
ATTN:  Library 

OEPARTHENT  OF  THE  NAVY 

Marine  Corps 

Department  of  the  Navy 
ATTN:  POM 

Marine  Corp  Dev  &  Education  Command 

Fire  Support  Section 

Department  of  the  Navy 

ATTN:  D091,  J.  Hartneady 

Naval  Construction  Battalion  Center 
ATTN:  Code  L08A 
ATTN:  Code  L51,  R.  Odell o 

Naval  Explosive  Ord  Disposal  Fac 

ATTN:  Code  504,  J.  Petrousky 

Naval  Postgraduate  School 

ATTN:  Code  1424  Library 

Naval  Research  Laboratory 
ATTN:  Code  2627 

Naval  Sea  Systems  Command 
ATTN:  SEA-9931G 
ATTN:  SEA-033 

Naya i  Surface  Weapons  Center 

White  Oak  Laboratory 

ATTN:  Code  U401,  M.  Kleinerman 
ATTN:  Code  F31 
ATTN:  Code  X211 

Naval  Surface  Weapons  Center 

ATTN:  Tech  Library  &  Info  Svcs  Branch 

Naval  Weapons  Center 

ATTN:  Code  266,  C.  Austin 
ATTN:  Code  233 

Naval  Weapons  Evaluation  Facility 
ATTN:  Code  10 

Office  of  the  Chief  of  Naval  Operations 
ATTN:  OP  S82 

ATTN:  OP  654C3,  R.  Piacesi 
ATTN:  OP  982E,  M.  Lenzini 

Strategic  Systems  Project  Office 

Department  of  the  Navy 
ATTN:  NSP-43 

DEI  An i fiENT  OF  THE  AIR  TQRCE 

Air  Force  Armament  Laboratory 

ATTN:  ADTC/XRS,  H.  Valentine 

3  cy  ATTN:  Dliv\  J.  Collins 

Air  Force  Institute  of  Technology 
ATTN:  Library 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 

Air  Force  Weapons  Laboratory 
Air  Force  Systems  Command 
ATTN:  SUL 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
ATTN:  INT 

Strategic  Air  Command/INT 
Department  of  the  Air  Force 
ATTN:  J.  McKinney 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
Department  of  the  Air  Force 
ATTN:  R.  Steere 

Foreign  Technology  Division 
Air  Force  Systems  Command 
ATTN:  NIIS  Library 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTM:  RSS 

Oklahoma  State  University 
Eld  Off  for  wpns  Effectiveness 
Department  of  the  Air  Force 
ATTN:  E.  Jackett 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN:  TSLD 

DEPARTMENT  OF  ENERGY 

Department  of  Energy 
Albuquerque  Operations  Office 
ATTN:  CTiD 

Department  of  Energy 
ATTN:  OMA/RD&T 

Department  of  Energy 
Nevada  Operations  Office 

ATTN:  Mail  &  Records  for  Tech  Lib 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Lawrence  Livermore  National  Lab 

ATTN:  Technical  Info  Dept  Library 
ATTN:  L-504,  M.  Wilkins 
ATTN:  J.  Goudreau 

los  Alamos  National  Scientific  Lab 
ATTN:  M/S632,  T.  Oowler 
ATTN:  MS  364 

Sandia  National  Laboratories 
Livermore  National  Laboratory 

ATTN:  Library  &  Security  Class  Div 


DEPARTMENT  Of  ENERGY  CONTRACTORS  (Continued) 
Sand i a  National  Lab 


ATTN: 

3141 

ATTN: 

W.  Caudle 

ATTN: 

A.  Chabai 

ATTN: 

J.  Colp 

ATTN: 

W.  Herrman 

ATTN: 

W.  Altsmeirer 

ATTN: 

W  Patterson 

ATTN: 

J.  Kaizur 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 

federal  Emergency  Management  Agency 
ATTN:  liazard  Eval  4  Vul  Red  Div 

NASA 

AMES  Research  Center 

ATTN:  R.  Jackson 

U.S.  Nuclear  Regulatory  Commission 

ATTN:  R.  Wliipp  for  Div  of  Sec  for  L.  Shao 

DEPARTMENT  0^  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  Technical  Information  Services 

Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc 
2  cy  ATTN:  0.  Trulio 

Avco  Research  &  Systems  Group 
ATTN:  D.  Henderson 
ATTN:  P.  Grady 
ATTN:  Library  A830 

BDM  Corp 

ATTN:  Corporate  Library 
ATTN:  T.  Neighbors 

Boeing  Co 

ATTN:  Aerospace  Library 

California  Research  &  Technology,  Inc 
ATTN:  K.  Kreyenhagen 
ATTN:  Library 

California  Research  4  Technology,  Inc 
ATTN:  D.  Orphal 

EG&G  Washington  Analytical  Services  Center,  Inc 
ATTN:  Library 

Engineering  Societies  Library 
ATTN:  A.  Mott 

General  Dynamics  Corp 
ATTN:  R.  Dibrell 

General  Electric  Company— TEMPO 
ATTN:  DAS1AC 

Honeywell,  Inc 

ATTN:  T.  Helvig 


DEPARTMENI  OF  DEFENSE  CONTRACTORS  fCnntim.pHl 

Institute  for  Defense  Analyses 
ATTN:  Classified  Library 

Kaman  AviDyne 

ATTN:  E.  Criscione 
ATTN:  N.  Hobbs 
ATTN:  Library 

Kaman  Sciences  Corp 
ATTN:  Library 

Lockheed  Missiles  4  Space  Co,  Inc 

ATTN:  Technical  Information  Center 

Lockheed  Missiles  4  Space  Co,  Inc 
ATTN:  It.  Culp 
ATTN:  TIC-Library 

Martin  Marietta  Corp 

ATTN:  H.  McQuaig 
ATTN:  A.  Cowan 
ATTN:  M.  Anthony 

Merritt  CASES,  Inc 
ATTN:  Library 
ATTN:  J.  Merritt 

Nathan  M.  Newmark  Consult  Eng  Svcs 
ATTN:  W.  Hall 
ATTN:  N.  Newmark 

University  of  New  Mexico 

Dept  of  Police  X  Parking  Security 
ATTN:  G.  Triandafal idis 

Pacific-Sierra  Research  Corp 
ATTN:  ii.  Brode 

Pacifica  Technology 
ATTN:  G.  Kent 
ATTN:  R.  Bjork 

Physics  International  Co 

ATTN:  Technical  Library 
ATTN:  L.  Behrmann 

R  &  D  Associates 

ATTN:  A.  Field 

ATTN:  Technical  Information  Center 

ATTN:  P.  Rausch 

ATTN:  W.  Wright,  Jr 

ATTN:  0.  Lewis 

ATTN:  P.  Haas 

Rand  Corp 

ATTN:  Library 

Science  Applications,  inc 

ATTN:  Technical  Library 

Science  Applications,  Inc 
ATTN:  W.  Layson 

SRI  International 

ATTN:  J.  Colton 
ATTN:  G.  Abrahamson 
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DEPARTMENT  Or  DEfENSE  CONTRACTORS  (Continued) 

Systems,  Science  S  Software,  Inc 
ATTN:  R.  Sedgewick 
ATTN:  Library 

Terra  Tek,  Inc 

ATTN:  Library 

TR1!  Defense  S.  Space  Sys  Group 

ATTN:-  Technical  Information  Center 
ATTN:  N.  Lipner 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Defense  &  Space  Sys  Group 
ATTN:  E.  Wong 
ATTN:  P.  Dai 

Weidlinger  Assoc,  Consulting  Engineers 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 

Weidlinger  Assoc,  Consulting  Engineers 
ATTN:  J.  Isenberg 


